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SUMMARY 


This  report  discusses  the  study  made  on  the  effectiveness  of 
various  matched  filters  for  processing  surface  waves  of  events  in  the  Kuriles 
as  recorded  at  LASA.  The  correlation  coefficients  and  signal-to-noise 
improvements  are  presented  for  matching  waveforms  which  include  master 
events,  a  chirp  waveform,  and  waveforms  generated  from  crustal  models. 
Based  on  the  results  of  this  study,  procedures  are  recommended  for  imple¬ 
mentation  of  large  scale  matched  filtering  at  the  SAAC  center. 
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SECTION  I 


INTRODUCTION 

The  purpose  of  this  study  was  to  investigate  the  practical  pro¬ 
blems  encountered  when  trying  to  implement  matched  filter  routinely  on  large 
numbers  of  events.  The  most  important  consideration  relative  to  large  scale 
implementation  is  the  choice  of  the  matching  waveform. 

In  choosing  a  matching  waveform,  the  net  gain,  in  terms  of  more 
useful  output,  must  be  balanced  in  terms  of  the  effort  required.  Any  procedure 
for  choosing  the  matching  waveform  must  be  logically  simple  enough  to  elimi¬ 
nate  the  need  for  human  judgement. 

The  answer  to  this  basic  question  is  considered  to  be  dependent 
upon  the  source  region  and  the  nature  of  the  crustal  path  between  source  and 
receiver.  This  study  was  conducted  using  17  earthquakes  in  the  Kuriles- 
Kamchitka  region  and  recorded  at  the  Montana  LASA.  Conclusions  drawn  from 
these  data  are  probably  not  applicable  to  other  source,  regions  and  recording 
stations;  however,  they  do  serve  as  an  initial  guide  to  wholesale  processing  of 
Rayleigh  and  Love  wavetrains. 

Two  basic  types  of  matching  waveforms  were  used  in  conducting 

this  study: 

•  Master  events  —  These  are  well  recorded  events 

from  the  same  source  region.  An  important 
consideration  is  how  close  the  maoter  event 
must  be  to  the  test  event,  as  this  governs  the 
storage  and  searching  required  for  general 
implementation.  The  master  events  were 
chosen  to  give  a  good  sampling  of  epicentral 
separations. 

•  Synthetic  waveforms  —  These  were  generated  from 

crustal  models.  A  simple  radar -type  Chirp 
(linear  increase  of  frequency  through  time  and 
uniform  amplitude)  was  also  used. 
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These  waveforms  were  then  used  to  do  matched  filtering  using 
the  ensemble  of  events.  Correlation  coefficients  and  signal-to-noise  (S/N) 
improvements  in  random  noise  were  calculated.  Finally  these  data  were  used 
to  make  recommendations  about  batch  processing  using  matched  filtering. 

The  information  contained  in  each  section  of  this  report  is: 

Section  II-brief  recapitulation  of  the  matched  filtering  concept;  Section  Ill- 
detailed  description  of  the  data  used  in  the  study  and  discussion  of  the  impli¬ 
cations  of  this  data  to  the  conclusions;  Section  IV-detail  discussion  of  the  re¬ 
sults  obtained.  The  conclusions  are  largely  subjective;  therefore,  the  experi¬ 
mental  results  are  presented  in  great  detail  for  the  interested  reader  to  examine; 
Section  V-conclusions  which  the  author  considers  pertinent  to  large  volume, 
automated,  matched  filtering  resulting  from  examination  of  the  data  and  results. 
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SECTION  II 


MATCHED  FILTERING  CONCEPTS 


Matched  filtering  is  a  convolution  of  a  matching  waveform  with 
the  signal-plus -noise  channel.  The  "best"  matching  waveform  is  the  time 
inverse  of  the  signal  waveform  itself.  It  can  be  shown  that,  under  the  most 
generally  accepted  criterion  of  goodness,  the  optimum  linear  filter  is 
shown  by  the  equation 


optimum  linear  filter  =  <!>  V 


(II- 1) 


where  4>  is  the  inverse  of  the  noise  correlation  matrix  and  V  is  the  signal 
waveform.  (This  assumes,  obviously,  that  a  prior  knowledge  of  the  signal 
waveform  exists). 

In  application  to  processing  of  seismic  surface  waves  the 
part  of  the  optimum  linear  filter  has  been  eliminated  and  some  estimate  of  the 
surface  wavetrain  has  been  used.  Studies  conducted  indicate  that  this  proced¬ 
ure  causes  a  negligable  loss  as  compared  with  the  optimum  linear  filter. 

For  random  noise,  which  has  the  same  bandwidth  as  the  signal, 
it  can  be  shown  that  matched  filtering  using  the  signal  waveform  which  is  of 
uniform  amplitude  would  give  a  gain  in  S/N  ratio  of  the  time  duration  -  band¬ 
width  product.  The  S/N  ratio  being  defined  as 

2 

S/N  _  (Max  signal  amp) _ 

Average  Mean  Square  Noise  (II-2) 


Lincoln  Lab  Test  Note  1967-50,  1967,  Long-period  signal  processing 
results  for  large  aperature  seismic  array,  Nov  15. 
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SECTION  III 

DESCRIPTION  OF  SUBJECT  EVENTS 


The  epicenter,  depth,  magnitude,  back  azimuth,  and  distance 
from  LASA  for  the  17  events  used  in  the  study  are  shown  in  Table  III  —  1 . 
Three  things  of  special  note  are 

•  All  the  events  are  near  magnitude  5.  This  should 
help  eliminate  great  differences  in  spectral  content 
due  to  source  characteristics,  and  should  give  a  good 
S/N  ratio  on  a  LASi  beam. 

•  All  the  events  are  shallow.  Deep  events  should  be 
eliminated  from  LP  processing  for  discrimination  by 
detection  of  depth  phases  from  short-period  stations. 

This  also  eliminates  another  possible  difference  in 
surface  wavetrain  characteristics  due  tc  source 
mechanisms  . 

•  All  the  events  have  a  back  azimuth  of  about  310  to 
315°  from  LASA.  This  means  the  surface  waves 

are  traveling  along  the  continent-ocean  basin  boundary 
over  most  of  their  path.  This  path  could  possibly 
result  in  anomalously  rapid  changes  in  waveform  as 
a  function  of  epicenter  (the  path  suggests  a  compli¬ 
cated  multipath  situation). 

The  geographical  locations  of  the  epicenters  of  the  events  are 
shown  in  Figure  III-l.  It  is  included  to  help  orient  the  reader  with  respect  to 
major  physiographic  features. 

The  LASA  beam  outputs  (vertical  traces)  for  each  event  are 
shown  in  Figure  III-2.  The  events  have  not  been  equally  scaled  in  this 
figure,  and  have  been  hand  traced.  This  imposed  some  unnatural  angularity 
to  the  waveforms. 
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Table  III- 1 


Event  summary 


Event 

Long 

}  Lat 

Depth 

Mb 

A 

Az 

1016 

150.  2E 

44.  6N 

33 

;  4.  7 

66.9 

311° 

1018 

140.  9E 

35.  6N 

59 

i  5-  3 

78.6 

310. 3° 

1019 

i  160. 4E 

52.  5N 

33 

4.  8 

56.  9 

312. 8° 

1131 

147.  2 

43.  5 

36 

5.7 

69.  2 

311.8 

1141 

147.  7 

43.  3 

12 

5.  3 

69.  1 

311.3 

1142 

147.  7 

43.  3 

33 

5.  1 

69.  1 

311.3 

2008 

152.  5 

46.  7 

40 

5.  6 

64.  3 

311.  5 

2010 

151. 7 

44.  3 

26 

5.  8 

66.  3 

310.  6 

2017 

151.4 

45.  5 

41 

5.  5 

64.  8 

311.4 

2018 

152.  0 

46.  3 

40 

5.  5 

64.  8 

311.4 

2019 

151.8 

45.  7 

40 

5.9 

65.  3 

311.8 

2020 

149.  5 

45.  2 

23 

5.4 

65.  3 

311.2 

2027 

155.  4 

47.  5 

29 

5.  2 

62.  3 

311.  1 

2028 

150.  7 

45.  4 

33 

5.  3 

66.  1 

311.4 

2029 

154.4 

49.  1 

33 

5.  3 

60.  7 

312.  7 

2032 

160. 4E 

54.  6N 

33 

5.  5 

55.  5 

315.  2° 

2033 

157.  7 

51.2 

24  f; 

5.  5 

58.  8 

313 

EVENT  2019 
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The  following  features  of  data  shown  in  Figure  III-2  are  con¬ 
sidered  pertinent: 

The  events  all  have  a  very  good  S/N  ratio. 

This  may  not  be  so  obvious  on  some  events 
as  the  time  tra  ;e  has  been  cut  off  and  shows 
only  a  gate  where  signal  energy  is  present. 

The  simplier  events  (such  as  2033  for  example) 
have  the  great  majority  of  their  energy  in  the 
first  500  seconds  after  onset.  Crustal  models 
suggest  that  at  60  to  70°  distance,  a  Rayleigh 
wavetrain  from  an  impulsive  source,  with  the 
bandwidths  observed  for  these  events,  should 
die  out  in  about  400  seconds. 

The  events'  waveforms  are  quite  variable. 

Many  of  the  events  are  highly  modulated  and 
many  have  a  considerable  coda  of  energy  which 
is  persistent  long  beyond  500  seconds.- 

Figure  III- 3  shows  the  amplitude  spectra  of  the  events  as  seen 
through  the  LASA  LP  instrumentation.  Shown  are  spectra  from  the  first  300, 
512,  and  1024  seconds  of  the  analysis  gate  which  indicate: 

•  With  the  exception  of  event  1019,  the  energy  in 
these  events  is  generally  limited  to  the  0.025- 
to  0.05-Hz  band  (20-  to  40-sec  periods). 

•  These  spectra  were  used  to  estimate  an  average 
bandwidth  for  each  gate  (300,  500,  and  1024  sec). 

This  average  number  listed  below,  was  used  for 
all  events  in  the  calculation  of  S/N  improvements 
contained  in  Section  IV: 

300  sec  - 0.015  Hz 

512  sec  - - 0.019  Hz 

1024  sec  — - 0.021  Hz 
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Figure  III-3.  A 
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Figure  III—  3 .  Amplitude  Spectra  of  Events  as  Seen 

Through  the  LASA  LP  Instrumentation.  (Sheet  1  of  2) 
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SECTION  IV 


RESULTS 

The  matched  filtering  was  done  using  events  1142,  2017,  2019, 
2029  and  2033  as  master  events.  These  events  were  chosen  to  give  a  good 
spread  of  epicentral  separations  when  used  as  the  test  event  with  each  of  the 
other  16  events.  Table  IV- 1  gives  the  epicentral  separation  between  the 
master  events  and  the  non-master  events.  The  separation  varies  from  0  to 
2216  km  and  a  fairly  good  sampling  of  epicentral  separations  was  obtained. 

Also  a  chirp  signal  was  used  as  a  matching  waveform  as  repre¬ 
sented  by  the  equation: 

Chirp  Signal  =  sin 

where: 


fo  + 


(fl  -  fQ> 
2T 


) 


(iv-1) 


fQ  =  0. 025  Hz 
fj  =  0. 05  Hz 

0  £  t  *  T  (Various  values  of  T  were  used) 

Finally,  surface  trains  derived  from  published*  dispersion 
curves  for  three  different  oceanic  models  were  used.  The  amplitude  spectra 
were  flat  from  0.  025  Hz  to  0.  05  Hz  and  tapered  off  to  zero  on  either  end  of 

this  band.  A  A  of  60  to  70°  was  used  for  each  model,  making  a  total  of  six 
such  Rayleigh  wavetrains. 

It  was  originally  intended  to  make  several  more  complicated 
models  using  continental  and  oceanic  in  various  reasonable  path  length  com¬ 
binations  and  assuming  no  phase  distortion  across  the  boundary.  The  results 


Saito,  M.  ,  H.  Takeuchi,  1966,  Surface  waves  across  the  pacific:  Bulletin 
of  the  Seismological  Society  of  America,  v.  56,  No.  5,  Oct.,  p.  1069-1078. 
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Table  IV- 1 


EPICENTRAL  DISTANCE  SEPARATION  (IN  KM)  BETWEEN 
MASTER  AND  OTHER  EVENTS 


1142 

2017 

2019 

2029 

2033 

1131 

44 

401 

439 

832 

1164 

1141 

0 

382 

421 

825 

1157 

1142 

0 

|  382 

421 

825 

1157 

2008 

531 

158 

123 

302 

628 

2010 

339 

135 

155 

572 

888 

2017 

380 

0 

38 

459 

786 

2018 

473 

100 

68 

359 

*  687 

2019 

418 

38 

0 

425 

750 

2020 

252 

152 

188 

570 

902 

2027 

759 

378 

340 

192 

444 

2028 

331 

55 

92 

497 

827 

2029 

822 

459 

425 

0 

331 

2032 

1154 

1197 

1163 

738 

419 

2033 

1154 

786 

750  ) 

331 

0 

1016 

244  j 

137 

175 

593 

922 

1018 

1062 

1437 

1475 

1884 

2216 

1019  | 

1388 

1018 

981 

567 

235 

using  the  above  models  were  so  discouraging  that  it  did  not  appear  likely  that 

any  one  crustal  model  would  work  appreciably  better  on  the  average  than  did 
Equation  IV- 1. 


Figure  IV- 1  shows  a  plot  of  the  correlation  coefficients  obtained 

using  the  first  300  seconds  of  both  the  master  event  and  the  test  event  using  the 
equation 


correlation  coefficients  = 


maxl  f  f?(t-+T) 

rp  I  1  A  1  C*  X  I 


(IV-2) 


The  correlation  coefficients  are  plotted  against  epicentral  separation  and  coded 
to  show  which  master  event  was  being  used.  Table  IV-2  lists  the  maximum 
value  of  each  of  the  correlation  functions  and  the  lag  at  which  it  occurred.  This 
same  data  is  listed  and  shown  for  512  seconds  in  Table  IV-3  and  Figure  IV-2 
and  for  1024  seconds  in  Table  IV-4  and  Figure  IV-3. 


Figure  IV-4  shows  the  correlation  coefficients  obtained  using  the 
various  synthetic  waveforms.  The  values  for  these  coefficients  are  listed  in 
Table  IV- 5. 
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Figure  IV- 1.  Correlation  Coefficients  of  Events  with  Master  Events  (First  300  Seconds). 
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Figure  IV- 2.  Correlation  Coefficients  of  Events  with  Master  Events  (First  512 
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Table  IV-6 

CORRELATION  COEFFICIENTS  AND  LAG 
USING  SYNTHETIC  WAVEFORMS  (FIRST  512  SECONDS) 
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The  features  of  this  data  are: 


•  The  correlation  functions  using  only  the  first 
part  (300  sec)  of  the  surface  wavetrain  are 
distinctly  higher  than  those  containing  later 
arriving  energy . 

•  A  simple  chirp  waveform  gives  a  correlation 
coefficient  almost  equal  to  those  obtained  using 
nearby  master  events  when  using  only  the  first 
300  seconds  of  the  events, 

•  The  correlation  coefficient  of  the  first  300 
seconds  of  the  Rayleigh  wavetrain  is  apparently 
stable  (chrip  works  well  and  master  events  work 
well  to  about  1000-km  epicentral  separation). 

The  correlation  coefficients  using  the  later  arriv¬ 
ing  energy  are  decidedly  lower,  and  degrade 
rapidly  with  epicentral  separation.  The  corre¬ 
lation  coefficients  using  the  entire  Rayleigh  wave- 
train  (1024  seconds)  are  low,  highly  scattered, 
and  indeed  some  of  them  are  not  meaningful  (i.  e.  , 
essentially  no  more  correlated  than  band  limited 
random  noise  could  reasonably  be).  This  last 
point  is  obvious  from  considering  the  lag  at  which 
they  occur  as  shown  in  Table  IV- 1. 

A  summary  of  the  values  of  the  average  correlation  coefficients 
given  in  Table  IV- 8  also  illustrates  these  tendencies. 


IV- 14 


•orvloos  group 


Table  IV -8 


AVERAGE  CORRELATION  COEFFICIENT 


Chirp  300  Sec 

0.  77 

Master  Waveform  300  sec  0-200  Km  Separation 

0.81 

Master  Waveform  300  sec  200-400  Km  Separation 

0.76 

Master  Waveform  300  sec  400-600  Km  Separation 

0.79 

Chirp  512  Sec 

0.648 

Master  Waveform  512  sec  0-200  Km  Separation 

0.76 

Master  Waveform  512  sec  200-400  Km  Separation 

0.  62 

Master  Waveform  512  sec  400-600  Km  Separation 

0.  59 

The  peak  of  the  matched  filter  output  can  be  used  directly  as  a 
discriminant  which  is  related  to  surface  wave  magnitude.  **  The  size  of  the 
peak,  for  a  fixed -amplitude  matched  filter,  is  dependent  upon  the  size  of  the 
arriving  surface  wave  and  the  correlation  between  the  wavetrain  and  the  match¬ 
ing  waveform.  For  the  discriminant  to  be  most  useful  it  should  be  related 
closely  to  the  size  of  the  arriving  surface  wavetrain.  For  estimating  surface 
wave  size  only  —  it  is,  therefore,  desirable  to  use  a  matching  waveform  which 
gives  as  high  a  correlation  coefficient  as  possible. 


Lincoln  Lab  Test  Note  1967-50,  1967,  Long-period  signal  processing  results 
for  large  aperature  seismic  array,  Nov.  15. 

** 

Seismic  Data  Laboratory  Report  No.  175,  1967,  Detection  of  surface  waves 
from  small  events  at  teleseismic  distances. 
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Based  on  the  events  studied, and  with  the  above  viewpoint  in  mind, 
the  following  matched  filtering  procedures  would  be  reasonable: 

•  Use  a  simple  chirp  waveform  of  length  somewhat 
less  (say  on  the  order  of  3/4)  than  the  expected 
length  of  a  Rayleigh  wave  train  from  an  impulsive 
source . 

•  Use  the  peak  of  the  matched  filter  output,  when 
statistically  significant,  directly  as  a  discrimi¬ 
nant  related  to  surface  wave  energy.  The  assump¬ 
tion  implicit  in  this  procedure  is  that  the  correlation 
coefficient  of  the  chirp  with  the  Rayleigh  v/avetrain 
is  the  same  for  each  event.  Based  on  the  events 
studied  this  does  not  seem  to  be  an  unreasonable 
assumption  for  the  first  300  seconds  of  the  data.  It 
may  be  that  the  above  assumption  would  not  be 
reasonable  in  comparing  events  from  widely  sepa¬ 
rated  locations. 

The  above  discussion  does  not  consider  the  S/N  improvement 
available  to  matched  filtering.  The  S/N  improvement  is  proportional  to  the 
time -bandwidth  product  of  the  signal.  It  may  therefore  be  a  good  strategy  to 
use  more  than  the  lirst  300  seconds  of  the  signal, which  means  accepting  a 
lower  and  more  variable  correlation  coefficient  but  gaining  improved  S/N 
ratio. 


The  S/N  ratio  was  defined  earlier  by  Equation  II-2,  and 

this  is  the  same  definition  used  in  previous  literature  dealing  with  matched 

filtering  of  seismic  signals.  Used  as  a  standard  for  comparison  of  S/N  ratio 

in  this  report  is  a  phaseless  bandpass  filter  which  passes  a  continuous  band 

containing  the  main  signal  energy,  and  removes  all  energy  outside  this  band. 

Such  a  filter  could  actually  be  approximated  very  simply  in  operation  by  using 

a  bank  of  bandpass  filters  and  taking  the  output  which  appears  to  give  the  best 
S/N  ratio. 


The  difference  between  the  bandpass  filter  and  the  matched 
filter  is  that  the  matched  filter  has  a  phase  which  is  designed  to  try  to  dupli¬ 
cate  that  phase  of  the  signal  which  causes  the  signal  energy  to  add  coherently 
and  peak.  Additionally,  like  the  zero-phase  filter,  it  is  a  bandpass  filter. 
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The  gain  of  the  matched  filter  over  the  zero-phaee  filter  is  a  function  of: 


•  Nonuniformity  of  the  wavetrain  (i.e.  ,  relationship 
between  peak-value-squared  and  mean-squared- 
value  ) 

•  Time  duration  of  the  signal  (if  only  the  first  part 
of  the  signal  is  used  some  potential  S/N  improve¬ 
ment  is  obviously  sacrificed) 

•  Bandwidth  of  the  signal 

•  Correlation  coefficient  between  the  matching  wave¬ 
form  and  the  subject  signal 

Figure  IV- 5  shows  the  S/N  improvements  relative  to  bandpass 
filtering  as  a  function  of  epicentral  separation  obtained  using  the  five  master 
events.  These  were  obtained  using  only  the  first  300  seconds  of  the  analysis 
■jate.  These  S/N  improvements  were  calculated  assuming  that  the  noise 
spectrum  was  flat  over  the  band  of  signal  energy  (about  0.  025  to  0.  05  Hz).  The 
numerical  values  of  the  numbers  plotted  in  Figure  IV-5  are  listed  in  Table  IV- 9. 

This  same  data  is  shown  in  figures  IV-6  and  IV-7  and  listed  in 
tables  IV -10  and  IV-11  for  the  512  and  1024  seconds  of  signals  respectively. 

Figure  IV-8  shows  the  S/N  improvements  obtained  using  a  repre¬ 
sentative  group  of  synthetic  waveforms  and  table  IV- 12  lists  the  numerical 
values.  Waveforms  which  were  inferior, or  not  appreciably  different  (on  the 
average), were  omitted  from  this  plot. 

The  S/N  improvements  that  could  have  been  obtained  with  the 
correct  matching  waveform  -the  signal  itself-  is  listed  in  Table  IV-1  3. 
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Table  IV- 9 


S/N  IMPROVEMENT  AS  A  FUNCTION  OF  EPICENTRAL  SEPARATION 

(FOR  30C  SECONDS) 


1142 

2017 

2019 

2029 

2033 

1131 

1.3 

1.0 

1.  1 

1.3 

1.3 

1141 

1.4 

1.  1 

1.2 

1.4 

1.4 

1142 

1.8 

0.  70 

0.  81 

1.2 

0.  81 

2008 

1.2 

0.  95 

1.  1 

1.3 

1.7 

2010 

0,94 

;  0.70 

0.78 

0.90 

0.91 

2017 

0.65 

1,7 

1.3 

1.2 

1.  1 

2018 

0.80 

1.2 

1.  1 

1.2 

1.2 

2019 

0.59 

0.96 

1.  3 

0.94 

0.  80 

2020 

0.88 

1.2 

1.  2 

1.2 

1.2 

2027 

0.69 

1.  1 

1.  2 

1.  1 

1.  1 

2028 

1.3 

0.  74 

0.  85 

1.2 

1.  1 

2029 

0.96 

1.0 

1.  1 

1.5 

1.  1 

2032 

0.  34 

0.  37 

0.25 

0.  35 

0.  61 

2033 

0.  74 

1.0 

0.97  ; 

1.2 

1.6 

1016 

0.  75 

1.2 

1.0 

1.  1 

1.  1 

1018 

0.61 

0.  17 

0.  25 

0.43 

0.44 

1019 

0.65 

1.  2 

1.5 

0.88 

0.67 
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SIGNAL- TO-NOISE  IMPROVEMENT  (POWER)  RELATIVE 
rO  WELL  CHOSEN  BANDPASS  FILTERING-LINEAR  SCALE 


SIGNAL -TO-NOISE  IM! 


Table  IV-  11 


S/N  IMPROVEMENTS  AS  A  FUNCTION  OF  EPICENTRAL  DISTANCES  USING 
MASTER  EVENTS'  SIGNALS  FOR  1024  SECONDS 


1142 

2017 

2019 

2029 

2033 

1131 

0.47 

1.  5 

2.  0 

1.4 

1.6 

1141 

0.  56 

0.  93 

1.8 

1.  1 

1.4 

1142 

4,  5 

0.  72 

0.  77 

0.63 

0.  85 

2008 

1.2 

1.4 

1.6 

1.3 

1.8 

2010 

0.29 

1.2 

1.  1 

0.74 

0.67 

2017 

0.98 

6.  1 

3.  1 

1.6 

1.6 

2018 

1.  1 

2.8 

3.9 

1.3 

0.93 

2019 

0.  53 

1.  55 

3.  1 

1.  15 

0.  81 

2020 

1.  1 

1.6 

3.4 

1.2 

1. 

2027 

0.86 

1.3 

1.5 

1.  1 

1.  : 

2028 

1.2 

1.3 

1.2 

0.  53 

1.6 

2029 

0.68 

1.3 

1.8 

4.9 

2.  2 

2032 

0.  2 

0.  32 

0.44 

0.29 

0.61 

2033 

0.  53 

0.73 

0.  73 

1.26 

2.  8 

1016 

0.  51 

1.  1 

1.7 

1.2 

0.61 

1018 

0.61 

0.44 

0.61 

0.72  | 

1.3 

1019 

0.  52 

1.4 

1.0 

0.42 

1.  1 
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S/N  IMPROVEMENTS  (POWER  NOT  IN  DB)  USING  THE  EVENT 
ITSELF  AS  THE  MATCHING  FILTER 


Event 

300  Seconds 

512  Seconds 

1024  Seconds 

1131 

1.7 

2.  6 

3.6 

1141 

1.9 

2.  5 

3.8 

1142 

j  1.8 

\  2. 8 

j  4.  5 

2008 

1.9 

3.4 

4.4 

2010 

1.4 

2.4 

3.  2 

2017q 

1.7 

3.7 

6.  1 

2018 

1  1.5 

3.  2 

4.9 

2019 

i.3 

2.  1 

3.  1 

2020 

1.6 

2.8 

4.0 

2027 

1.5 

3.  1 

3.6 

2028 

1.8 

2.9 

4.  1 

2029 

1.  5 

3.  5 

4.  9 

2032 

1.3 

1.8 

2.  1 

2033  ; 

1.6 

2.  5 

2.8 

1016 

1.6 

rsj 

• 

oo 

5.  1  j 

1018 

1.7 

2.  8 

3.4 

1019 

1.7 

3.6 

4.7 

Uniform  Amplitude 
Event  With  Same 
Bandwidth 

4.4 

9.6  { 

21.2 

These  data  have  several  important  features: 

*  The  signals  do  not  offer  the  possibility  of  large 
gains  m  S/N  ratio.  It  can  be  seen,  by  the  data 
listed  in  table  IV- 13,  that  the  S/N  improvements 
available  with  perfect  matching  average:  1.6  (2db) 
for  300  seconds,  2.  85  (4.  5db)  for  512  seconds, 
and  4.0  (6db)  for  1024  seconds.  The  1024  seconds 
includes  essentially  all  the  signal  energy  in  every 
case.  These  low  figures  occur  basically  because 
the  square  of  the  peak-signal-excursion  is  large 
relative  tc  twice  the  mean-square-signal  value. 

For  a  uniform  wavetrain  the  square  of  the  maxi¬ 
mum  excursion  would  be  about  twice  the  mean 
squared  value  of  the  wavetrain.  This  would  re¬ 
sult  in  a  S/N  improvement  of  about:  4.  4  (6.  4db) 
for  the  first  300  seconds  of  data  or  9.  6  (9.  8db)  for 
512  seconds  of  data 

•  The  S/N  results  actually  realized  using  the  first 
300  seconds  of  data  are  not  much  different  than 
bandpass  filtering  on  the  average 

•  In  order  to  use  some  of  the  later  arriving  energy 
(512  or  1024  seconds),  it  is  necessary  to  have  a 
master  event  very  near.  Otherwise  there  is  a  good 
chance  of  considerable  S/N  loss.  For  this  data 
the  average  S/N  improvement  factor,  using  master 
events  within  200  km,  was:  1.67  for  the  first  512 
seconds  and  1.71  for  the  first  1024  seconds.  There 
were  no  examples  of  S/N  loss  using  512  points  and 
only  2  examples  of  loss  using  1024  points.  The 
synthetic  waveforms  give  significantly  poorer  re¬ 
sults  than  nearby  master  events  on  the  512-second 
and  especially  on  the  1024-second  data 

•  Unlike  the  data  using  only  300  points,  the  S/N  im¬ 
provements,  using  master  events  and  a  longer  piece 
of  the  wavetrain,  decay  rapidly  with  increase  in  the 
epicentr al  separation.  It  appears  to  be  very  difficult 
to  obtain  the  S/N  improvements  inherently  available 
in  a  elongated  wavetrain.  (e.  g.  ,  A  6.  1(7.  9db)  gain 
would  be  possible  with  an  exact  match  to  the  entire  wave- 
train  of  event  2017,  but  the  master  events  give  the 
improvements  shown  in  Table  IV- 14. 
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Table  IV- 14 


MASTER  EVENTS  S/N  IMPROVEMENTS 
FOR  EVENT  2017 


Master  Event  No. 

1142 

2019 

2029 

2033 


Epicentral  Separation 

Gain  Factor 

(km) 

(Linear) 

380 

0.98 

38 

3.  1 

459 

1.6 

786 

1.6 
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CONCLUSIONS 


Although  redundant,  it  is  felt  that  it  is  worthwhile  to  repeat  the 
significant  features  of  the  signals  studied  before  presenting  the  conclusions. 
These  features  are: 

•  Events  were  recorded  at  LASA  from  the  Kuriles 
region  and,  therefore,  traveled  a  path  along  the 
continent-ocean  basin  boundary.  This  path  might 
well  be  subject  to  complex  multipathing,  thereby 
resulting  in  unusually  modulated  (nonuniform) 

wavetrains  and  wavetrains  which  varied  quite 
rapidly  as  a  function  of  epicenter.  Both  of  these 
effects  are  seriously  detrimental  to  matched  filtering. 

To  illustrate  this  further  group  velocity  curves  cal¬ 
culated  from  events  2017  and  2018  are  shown  in 
Figure  V-l. 

Not  only  are  these  group  velocity  plots  complex  but 
they  are  rather  different,  even  though  the  events 
are  separated  by  only  100  Km.  Many  of  the  events 
gave  a  group  velocity  plots  which  were  more  scattered 
and  complex  than  these. 

By  contrast,  a  group  velocity  plot  is  shown  for  an 
event  from  the  South  Pacific.  For  events  from 
this  region  and  most  other  regions,  the  group 
velocity  plots  cluster  along  a  smooth,  respectable 
curve. 

•  The  events  studied  were  generally  band  limited  to  the 
0.  025-  to  0.05-Hz  range.  However,  other  studies* 
indicate  that  LP  surface  phases  routinely  recorded 
have  a  somewhat  broader  bandwidth.  In  such  a  case 
the  increased  S/N  improvement  available  to  match 
filtering  is  directly  proportional  to  the  signal  bandwidth. 

Pomeroy,  et  al,  1969,  Preliminary  results  from  high-gain,  wide-band, 
long-period  electromagnetic  seismograph  systems:  Journal  of  Geophysical 
Research,  v.  74,  No.  12,  June  15. 
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The  following  conclusions  are  made  as  a  result  of  data  presented  in  this  report. 

For  the  signals  studied  from  the  Kuriles,  routine  matched  filter¬ 
ing  at  LASA  does  not  appear  to  be  a  very  powerful  tool.  Two  routine  methods 
of  operation  appear  reasonable. 

Since  the  output  of  the  matched  filtering  operation  can  be  used 
directly  as  a  discriminant  related  to  the  total  energy  in  the  wavetrain’  (in  a 
fixed  region  at  least),  it  would  be  desirable  to  make  this  measurement  using  a 
matching  filter  which  vould  give  a  high  correlation  coefficient  with  most  events. 
The  signals  analyzed  indicate  that  this  could  be  accomplished  by  using  only 
aoout  3/4  of  the  expected  duration  of  a  simple  event  and  employing  either  a 
.<airly  close  (800  km)  master  ever.t  or  a  simple  chirp  waveform.  Both  of 
these  gave  average  correlation  coefficients  on  the  order  of  0.8.  There  were 
not  sufficient  signals  processed  to  give  a  good  estimate  of  the  variability  of 
such  a  procedure,  but  an  assumption  of  equal  correlation  with  each  event  does 
not  appear  to  be  unreasonable. 

This  procedure  would  only  give  S/N  improvement  comparable 
to  a  narrow  band  filter,  but  it  would  give  an  output  which  would  be  more 
useful  for  discrimination  purposes.  It  will  be  recalled  that  the  bandpass  filter 
does  not  compress  the  waveform;  therefore,  a  comparable  S/N  ratio  means 
that  the  maximum  excursion  of  the  bandpassed  signal  would  be  about  as  visible 
in  the  noise  as  the  peak  of  the  matched  filter  output.  A  characterization  of  the 
surface  wavetrain,  using  the  matched  filtered  peak,  would  generally  be  preferred 
to  information  on  the  size  of  the  maximum  excursion. 


Seismic  Data  Laboratory  Report  No.  175,  1967,  Detection  of  surface  waves 
from  small  events  at  telesuismic  distances. 
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If  it  is  desired  to  improve  upon  the  above  procedure  it  would 
require  a  close  (within  about  200  km)  master  event.  Using  such  a  master 
event  and  the  entire,  or  nearly  the  entire,  wavetrain,  a  small  S/N  improve¬ 
ment  (a  factor  of  about  1.7  or  2.  3  db)  was  obtained.  The  average  correlation 
coefficient  was  0.755  for  512  seconds  using  master  events  within  200  km.  This 
is  comparable  to  what  could  be  obtained  using  a  chirp  wave  form  and  300  seconds 
(0.  77).  Again,  there  was  insufficient  data  to  make  a  really  meaningful  estimate 
of  the  variability  of  these  correlation  coefficients. 

It  is  recommended  that  the  above  procedures  be  implemented 
initially  in  the  off-line,  long-period,  analysis  package  at  SAAC.  After  an 
accumulation  of  a  more  significant  data  base  these  procedures  should  be  re¬ 
evaluated  on  a  regional  basis. 
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